Origin of synchronized trafRc flow on highways and its dynamic phase transitions 
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We study the traffic flow on a fiighway witli ramps tlirough numerical simulations of a hydrody- 
namic traffic flow model. It is found tliat tfie presence of the external vehicle flux through ramps 
generates a new state of recurring humps (RH). This novel dynamic state is characterized by tem- 
poral oscillations of the vehicle density and velocity which are localized near ramps, and found to be 
the origin of the synchronized traffic flow reported recently [PRL 79, 4030 (1997)]. We also argue 
that the dynamic phase transitions between the free flow and the RH state can be interpreted as a 
subcritical Hopf bifurcation. 
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Experiences show that traffic flow has complicated 
properties. The fact that automobile is one of main trans- 
portation tools raises the trafRc flow as one of the most 
important problems for engineers 0|. For physicists, on 
the other hand, traffic flow is an interesting many-body 
problem of interacting vehicles. Numerous experimental 
measurements revealed that trafRc flow possesses quali- 
tatively distinct dynamic states In particular, three 
distinct dynamic phases are observed in highways : the 
free trafRc flow which is analogous to laminar flow in 
fluid systems, the traffic jam state where vehicles almost 
do not move, and the synchronized traffic flow which is 
characterized by complicated temporal variations of the 
vehicle density and velocity. 

Paralleled with experiments, many physical models 
have been proposed Cellular automaton models H] 
have been developed which simulate each individual vehi- 
cle, and hydrodynamic models which provide macro- 
scopic description of traffic flow. Subsequent studies (^H] 
of the models have explained many observed features of 
the free flow and traffic jams in highways. However, 
no satisfactory explanation for the synchronized flow is 
available to our knowledge. 

Recently, Kerner and Rehborn reported analysis of 
systematic measurements performed on German high- 
ways. As one of main results, it was pointed out that 
the synchronized flow is spatially localized near ramps 
on highways This observation motivated us to ex- 
plore in this paper effects of ramps on highway traffic 
flow. Through numerical simulations of a hydrodynamic 
model, we find that the presence of ramps generates a 
new kind of traffic states which becomes a spatially local- 
ized limit cycle of highway traffic flow under the constant 
external flux. We examine properties of the novel state 
and show that it is the origin of the synchronized flow. 

In this work, we adopt the hydrodynamic model of 
highway traffic flow proposed by Kerner and Konhauser 
, where the dynamic evolution is described by the 



Navier-Stokes-type equation of motion, 




Here p{x, t) is the local vehicle density, v{x, t) the lo- 
cal velocity, V{p) the safe velocity that is achieved in 
the time-independent and homogeneous traffic flow, and 
T,co,fi are appropriate constants. Eq. is paired with 
the modified equation of continuity 



'dt dx ^ Qin{t)f{x - Xin) - qout{t)ip{x - Xout) , 

(2) 

where the source and the drain terms on the right hand 
side represent the external flux through an on-ramp and 
through an off-ramp, respectively p^ . Here 1^9(2;), de- 
scribing the spatial distribution of the external flux, is lo- 
calized near x — and normalized so that qin{t) iloutit)) 
represents the total incoming (outgoing) flux. 

To study effects of a single ramp, two ramps are 
separated by a large distance {\xm — Xout \ = L/2 where 
L is the system size) and numerical simulations are per- 
formed with periodic boundary conditions. The two-step 
Lax-Wendroff scheme is adopted as the main simulation 
scheme, and its reliability is verifled through comparison 
with an alternative scheme: the classical fourth-order 
Runge-Kutta scheme applied to the time and the cen- 
tered Euler scheme to the space [Q. Simulations are 
carried out for many different sets of parameters and 
qualitatively the same results are obtained. So for defl- 
niteness, we present results only for the following choice 
of parameters : t ~ 0.5 min, /i = 600 km/h, cq = 54 
km/h, and V{p) = Voil - p/p)/(l + E{p/pY) where 
the maximum density p — 140 vehicles/km, Vq = 120 
km/h, E = 100, and — A ||l^. Concerning the dis- 
cretization, spatial intervals of Ax — 37.8 m and time 
intervals of A< = 10^'' min are found to be suitable. 
We choose the spatial distribution of the external flux as 
ip{x) = (27rcr2)-i/2 exp {-x"^ /2a'^) with a = 56.7 m. 
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This model was investigated previously for 

the constant external flux qi„{t) = qout{t) = /■ For 
small /, it was found that an initially homogeneous flow 
p{x,0) = ph, v{x,0) = V{ph) evolves to slightly modi- 
fied free flow where homogeneous regions with different 
densities are separated by narrow density-rising (or de- 
scending) regions near the ramps, the so-called transition 
layers. In contrast, for / larger than a critical value, a 
local avalanche- like process occurs at the transition layer 
and traffic jam appears spontaneously. This study, how- 
ever, failed to probe the synchronized traffic flow. 

To find a clue to the missing third phase in traffic flow, 
we pay attention to the experimental observation ||ic| ] 
that for a range of /, traffic fiow can be either in the 
synchronized flow or in the free flow. This bi-stability 
suggests that the transition from one locally stable state 
to the other may require some triggering events. So in 
our simulations, we apply a pulse-type perturbation with 
a finite amplitude. Specifically we first prepare a transi- 
tion layer by applying the constant external fiux / which 
is below the critical point fc (for f > fc, the stable free 
fiow does not exist) . Then a pulse of additional flux Sqin 
is applied at the on-ramp for a short duration St. As a 
result, a localized oscillating state appears from the free 
fiow @ (Fig. 0). 
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FIG. 1. Birth and evolution of the RH state from the free 
flow is shown in this plot of the density proflle (p — pupstream) 
near the on-ramp (xi„=18.9 km) for the average density of the 
system ph=22.4 vehicles/km (?s {pupatream + pdo wnstream )/2) 
and the constant external flux / = 318 vehicles/h (the sys- 
tem size L = 75.6 km). Transition from the free flow to the 
RH state is triggered by a pulse-type perturbation applied at 
t — 50 min with (5gi„=318 vehicles/h, St — 5 min. 



After a transient period, the localized oscillation be- 
comes periodic in time. We observe that the Fourier 
spectrum of the oscillation shows sharp peaks at each in- 
teger multiple of the basic frequency 1/T (« 0.068 min^^ 
in case of Fig. 0). In the simulation, it turns out that 
properties of this periodic asymptotic state, such as the 
period and the oscillation amplitude, are essentially in- 
dependent of Sqin and St as long as they are large enough 
to trigger the transition. This strongly suggests that the 



periodic oscillation is not a transient process with long 
decay time but a limit cycle of Eqs. (0,§). For definite- 
ness, we call this state of traffic flow the recurring hump 
(RH) state in this paper, which will be later compared 
with the synchronized flow. 

Limit cycles are generated in numerous examples of 
nonlinear autonomous systems (that is, systems without 
explicitly given time dependence) |17jl^. In the present 
example of traffic fiow, the limit cycle can be character- 
ized as a self-excited (autocatalytic) oscillator (see Sec. 
5.6 in Ref. [p^), where constant external fiux serves as 
a source of periodically generated excitations (humps). 
Excitations are, however, relaxed within a localized re- 
gion. When the upstream vehicle density is lower than 
the critical value Pcr{~ 25 vehicles/km for our parame- 
ter choice) , a localized inhomogeneity decays away in the 
homogeneous traffic environment unless the amplitude of 
the inhomogeneity is larger than a critical magnitude [|| . 
Thus humps cannot survive far away from the on-ramp 
if its size does not exceed the critical magnitude. In this 
way, the localization can be achieved. 

The character of the localized oscillation becomes ev- 
ident in the density-flow diagrams {p{x,t) vs.q{x,t) = 
p{x,t)v{x,t)). In contrast to a straight line for the free 
fiow, the density-flow relation for the RH state forms a 
closed loop at a; = Xin (Fig. ||(a)), which implies the peri- 
odicity of the oscillation and also the phase difference in 
oscillation between p(x, t) and ^(x, t). As a: moves down- 
stream, the loop deforms gradually to a smaller loop and 
eventually joins the free flow (Fig. |^(a)), which is a con- 
sequence of the localization. We also examine the effect 
of randomly fluctuating external flux. Figure |^(b) shows 
that though the exact periodicity is lost, the oscillation 
itself is still stable under random fluctuations. 

Below we investigate the transition from the free flow 
to the RH state. For definiteness, we fix the perturba- 
tion, Sqin = 159 vehicles/h, St = Q min and apply it to 
the transition layer generated by / (Fig. ^(a)). For small 
/, the free flow survives the perturbation. For / larger 
than a critical value /i, however, the /jnite-amplitude RH 
state is induced. We emphasize that a finite perturba- 
tion is essential for the transition. As the perturbation 
becomes weaker, /i becomes larger and for SqinSt —>■ 0, 
the transition to the RH state does not occur for the 
whole range of / smaller than fc [|ll|. For the backward 
transition, on the other hand, it turns out that it may 
occur even without flnite amplitude perturbations. As 
/ decreases adiabatically from / > /i, the amplitude of 
the RH state varies as in Figure ||(a). The system first 
follows its old path. Below / = /i, however, the system 
still remains in the RH state instead of going back to the 
free fiow. The RH state is maintained until / reaches a 
lower critical value /2, where the transition to the free 
flow occurs We mention that transitions between 

the RH state and the free flow show the same hysteresis 
as measured in highways pO[ ] . The influence of the transi- 
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tion on traffic flow becomes clear in the following natural 
order parameter: the spatio-temporal average velocity, 



(V) EE — 



t+T 



dt' 



R./2 
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dx v{x + Xin,t') , 



(3) 



where T is the period of the RH state and R is the size 
of the averaging range. In Figure |^(b), (u) makes discon- 
tinuous jumps at the transition points /i and /2 |20|]. 
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FIG. 2. (a) Density-flow diagrams for the average den- 
sity ph = 22.4 vehicles/km, / = 318 vehicles/h measured 
at 2: — Xin = km (the largest loop), 0.9 km, and 3.8 km (the 
smallest loop), respectively. The dashed line represents the 
Iree flow diagram, (b) Density-flow relation at the on-ramp 
under the presence of random fluctuations in the external flux. 
Data are shown for qin.out{t) = 318 vehicles/h +5qin,out{t) 
where Sqin,out{t) is a random constant function in each time 
interval of Tp = 5 min. At the end of each interval, the values 
of two independent random fluctuations 5qin{t),5qout{t) are 
reset by new random numbers uniformly distributed in the 
range [—95,95] vehicles/h. 

Another interesting property of the RH state appears 
in multi-lane situations. To demonstrate this, we ex- 
tended the traffic equations to a two-lane system. The 
equation of motion (0) and the continuity equation (|^) 
apply to each lane i — 1,2. We assume that ramps are 
connected to the lane 2 and so the source and drain terms 
appear in the continuity equation for the lane 2 only. 
We simulate the inter-lane interaction effect in a mini- 
mal way by introducing to the continuity equations lane- 
change terms {dpi/ dt)ch that account for the inter-lane 
flux due to the lane change of vehicles. For the simple 
choice, {dpi/dt)ch = a{pj{x,t) - pi{x,t)) {i ^ j), it is 
found that when the flow in the lane 2 makes the transi- 
tion to the RH state, it is accompanied by appearance of 
the synchronized oscillations of the velocity and the den- 
sity in the lane 1 (Fig. ||(c)). This property of the syn- 
chronization is examined for different functional forms of 
{dpi/dt)ch as well since its precise form is not well de- 
termined yet. Qualitatively same results are recovered 
in all cases, which demonstrates that the synchroniza- 
tion is a generic property of the RH state in multi-lane 
situations |21|. It is worth commenting that this kind 
of synchronization phenomena is a common property in 
many examples of self-exciting systems (see Sec. 5.13 in 
Ref. P). 



Now it should be noticed that many properties of the 
RH state are identical to those of the synchronized flow 
jlOf , e.g., the discontinuous transition from the free flow 
to the synchronized flow induced by localized perturba- 
tions of finite amplitudes, hysteresis, stability of synchro- 
nized flow (hours of self-maintenance), gradual spatial 
transitions from synchronized flow to free flow, and syn- 
chronized oscillations. Therefore, we conclude that the 
RH state is the origin of the synchronized traffic flow. 
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FIG. 3. (a) Amplitude of density oscillation measured at 
the on-ramp for ph = 22.4 vehicles/km. The zero amplitude 
implies free flow. Discontinuous jumps occur at /i ~ 241 
vehicles/h and /2 ^ 184 vehicles/h. (b) {v) (Eq. ^) as a 
function of / with R = 7.6 km. Though the precise value of 
{«) depends on R, the presence of discontinuities is universal, 
(c) Synchronized oscillations in a two-lane system. The figure 
shows the temporal variations of the velocities at the on-ramp 
for the lane 1 (solid line) and for the lane 2 (dotted line). A 
pulse of the external fiux introduced at t = 50 min induced 
the transition to the RH state. Notice that the variations 
are synchronized in both lanes, (d) Schematic diagram of the 
bifurcation scenario. The point A corresponds to a turning 
point and B to a subcritical Hopf bifurcation point. While the 
transition at /2 is spontaneous, the transition at /i should be 
aided by an external triggering event. As a result, the value of 
/i depends on the strength of triggering events. For / > /c, 
the free flow loses its stability spontaneously. 



We interpret our results within the standard frame- 
work of nonlinear dynamics. The free flow corresponds 
to a point attractor. On the other hand, many features of 
the RH state, such as stability and discontinuous tran- 
sitions, assert that the RH state corresponds to a sta- 
ble limit cycle. Hysteresis implies bi-stability for certain 
range of /. The spontaneous backward transition from 
the limit cycle to the point attractor means that the lower 
end of the bi-stable region is f^- On the other hand, the 
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necessity of finite perturbations for the forward transi- 
tion suggests that the upper end of the bi-stable region 
goes above /i and in fact extends to fc- Discontinuous 
transitions in both directions indicate that the limit cy- 
cle is not connected to the point attractor (Fig. ||(c)). 
The discontinuity also suggests that there should exist 
still another asymptotic state, which serves as a bound- 
ary between the basins of the attraction toward the two 
locally stable asymptotic states. One plausible candidate 
for the boundary is an unstable limit cycle that is con- 
nected to the stable limit cycle at / = /2 and also to 
the fixed point at f — fc (Fig. ||(c)). In this case, the 
whole transition behaviors are results of a turning point 
(/2) combined with a subcritical Hopf bifurcation (fc)- 

Lastly we briefly discuss the synchronized flow far away 
from ramps reported in Refs. There is one im- 

portant difference between the synchronized flow near 
ramps and that far away from ramps. While the former 
is a non-decaying state of oscillations, the latter appears 
only as a transient process; after some upstream move- 
ment since its creation, it either disappears or transforms 
into a jam pO[ |. Therefore, the synchronized flow far away 
from ramps is not a stable dynamic phase of traffic flow, 
and is not studied in this paper. The appearance of this 
transient process requires further investigations. 

In summary, we find that there exists recurring hump 
(RH) state in highway traffic flow with ramps. In this 
state, the density and the flow oscillate periodically and 
the oscillations are localized near the on-ramp. The RH 
state is a stable limit cycle of the nonlinear traffic equa- 
tions (|l|,||). The transition between the free flow and the 
RH state is discontinuous and shows hysteresis. Many 
features of the RH state are identical to those of the syn- 
chronized flow and thus we conclude that the RH state 
is the origin of the synchronized flow observed in real 
highways 
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